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ABSTRACT 

The most massive galaxies in the Universe are also the oldest. To overturn this apparent contra- 
diction with hierarchical growth models, we focus on the group-scale haloes which host most of these 
galaxies. Our z ~ 0.4 group sample is selected in redshift space from the CNOC2 redshift survey. A 
stellar mass selected M* > 2 x 1O 1O M0 sample is constructed using IRAC observations. A sensitive 
Mid InfraRed (MIR) IRAC colour is used to isolate passive galaxies. It produces a bimodal distribu- 
tion, in which passive galaxies (highlighted by morphological early-types) define a tight MIR colour 
sequence (Infrared Passive Sequence, IPS). This is due to stellar atmospheric emission from old stellar 
populations. Significantly offset from the IPS are galaxies where reemission by dust boosts emission 
at A bs=8/im. We term them InfraRed-Excess galaxies whether star formation and/or AGN activity 
are present. They include all known morphological late-types. Comparison with EW[OII] shows that 
MIR colour is highly sensitive to low levels of activity, and allows us to separate dusty-active from 
passive galaxies at high stellar mass. The fraction of InfraRed Excess galaxies, f(IRE) drops with 
M*, such that f(IRE) = 0.5 at a "crossover mass" of M cr ~ 1.3 x 10 n M Q . Within our optically- 
defined group sample there is a strong and consistent deficit in f(IRE) at all masses, but most clearly 
at M* > 10 x1 Mq. Suppression of star formation must mainly occur in groups. In particular, the 
observed trend of f(IRE) with M» can be explained if suppression of M* > 10 11 M Q galaxies occurs 
primarily in the group environment. This is confirmed using a mock galaxy catalogue derived from 
the Millenium Simulation. In this way, the mass-dependent evolution in f(IRE) {downsizing) can 
be driven solely by structure growth in the Universe, as more galaxies are accreted into group-sized 
haloes with cosmic time. 

Subject headings: galaxies: statistics — galaxies: evolution — infrared: galaxies — galaxies: photom- 
etry — galaxies: high-redshift — galaxies: clusters: general 



1. INTRODUCTION 

Perhaps the biggest puzzle of galaxy formation in 
recent decades has been the early and rapid forma- 
tion of the most m assive galaxies ()Bender et al.l 119981 : 
iCimatti et al.1 l2004h . This apparently suppo rted the 
"monolithic collapse" model (|Eggen et al.lll962T ). in con- 
tradiction to the standard picture of a Cold Dark Matter 
(CDM) Universe, in which galaxies grow hierarchically 
through mergers, along with t he dark matter hal oes in 
which they are embedded (e.g. iBenson et ai1l2002f ). Ob- 
servationally, the picture has been embellished by re- 
cent galaxy surveys from z ~ to z ~ 1.5. Going to 
higher galax y mass and to denser environments, galax- 
ies are older (iKauffmann et al.ll200l lPanter et all 1 120031: 
Poggianti et al.l 120041 : iThomas et al l l2005t ISalim et al l 



2005) and more likely to have exper ienced a shorter star 
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downsizing ( Cowie et al. 


11999ft is used to describe the 



earlier and less extended star formation histories in more 
massive galaxies. 

Recent attempts to explain downsizing trends have ex- 
amined the ways in which gas heating and cooling depend 



upon mass and redshift (e.g. iDekel fc Birnboid 1200(1 . 
and invoked new heating methods such as feedback 
from Active Galac tic Nuclei (AGN. [Bower et al.1 [20061 : 
ICroton et al.ll2006l ). Regardless of the physics, a hierar- 
chical Universe can itself supply a natural explanation 
for dependence of galaxy properties up on stellar mass 
and u pon environment at the same time (|De Lucia et al.l 
2006). As more and more massive galaxies are encor- 
porated into more and more massive haloes, the prob- 
ability of star formation being truncated is increased. 
Even in the populous loose group environment the prop- 
erties of galaxies depend st rongly upon environment (e.g . 
Postman fc Gelid Tl98l IZ abludof f fc Mulchaevl 119981 ; 
Lewis et a l. 2002; Gom ez et al.ll2 003: Bal ogh et aflliZOG 7 
Wilman et al.l 120051 : IWeinmann et al.ll2006t iGerke etal 
20071) . Thus, evolution in groups can play an important 



role in driving evolutionary trends. 

The fraction of passive galaxies increases steadily and 
monotonically with increasing stellar mass. It is there- 
fore convenient to quantify downsizing trends by iden- 
tifying a stellar mass at which the fraction of passive 
(or red/early- type) galaxies is defined to be exactly 
50%. This has no physical significance, and so we pre- 
fer to call this the "crossover mass" (M cr ) rather than 
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the 11 transition mass " (e.g. iBundv et all I2006H . M cr 
decr eases strongly with increasing environmental den- 
sity (iBaldrv et al"j|2006t IBundv et aTll2006l; lHaines et all 



2007), with a maximum M cr - 2 x 10 1U M in the lowest 
density environments wher e the z ~ red galaxy fraction 
is 50% (|Baldrv et alJl2006f >. The evolution of M cr for the 
global population, indicates M cr ~ 2.5 — 6 x 10 10 M 
at z ~ 0.4 and M cr - 7 - 10 x 10 10 M at z ~ 1 
(|Hopkins et al.ll2007ft . 

These studies, mostly at rest-frame UV to optical 
wavelengths, are hindered by the strong correlation be- 
tween gala xy mass and dust attenu at ion in star-forming 
galax ie s dGiovanell et all 119951: I Wang fc Heckmanl 
19961: iMasters et all \2Q03L iBrinchmann et all 12004 



Weiner et al.ll2007[ ). This correlation may arise directly 



from the variations of star formation history with mass 
(|Calura et al.l I2007D and can lead to confusion between 
truly passive and dusty star forming galaxies. In some 
cases light from the youngest stars can be totally re- 
moved from the line of sight, not only because their SED 
peaks at shorter wavel engths where dust attenuation 
is most effective (e.g. iPierini et al.l I2004D . but also 
because the HII regions i n which they ar e embedded 
can b e particularly dusty (|Duc et alj l2002t I Tuffs et all 
2004). This degeneracy of passive versus dusty star 
forming galaxies is well known at high redshifts, where 
a selection of optically red galaxies (Extremely Red 
Object/Galaxy -ERO/ERG or Distant Red Galaxy - 
DRG) is more and more likely to select dusty starbursts 
as one moves to even higher redshift (from ~ 50% at 
z ~ 1 (ICimatti et ail 120021: lLotz et alJ [20Qgft to > 90% 
at z ~ 2-3 (jPanovich et al.ll2006h . with an increasing 
contributi on from increasingly active infrar ed-bright 
galaxies ([Caputi et alJl2006t iDaddi et al.ll2007h ). 

Even in th e local Universe the optical red sequence (e.g. 
iBower et al.l [1992) is contaminated with star-forming 
galaxies and some AGN. These "interlopers" can be 
identified using optical or UV colour-colour diagrams 
or high signal to noise emission line measurements (e.g. 
I Wolf et alJfeOOa lHaines et al.ll2007l ). However the dust 
emission really stands out in the infrared where the ab- 
sorbed e nergy (mainly non-i onizing UV photons) is re- 
emitted (jPopescu et al.1 12000L and references therein) . 
31% of red-sequence galaxies (by number density) are 
24/im-bright, corresponding to ~ 17 % star- forming 
galax ies and ~ 14% AGN-host galaxies ([Davoodi et al.l 
l200l . 

At shorter wavelengths (^3-12/^m) strong emis- 
sion features exist in t h e presence of sta r forma- 
tion ([Phillips et al.l 11984 iRoche et al] Il991h . How- 
ever they are usually absent in th e more highly ion- 
ized environments typica l of AGN (jAitken et alj [19821 
IRoche et all 11984 I1991D . These have been identi- 
fied as PAH (Polycyclic- Aromatic- Hydrocarbon) features 
(most obviously characterized by the strong, broad 
emission features at 6 .2/zm 7.7/xm 8.6Atm and 11.3/im, 
iLeger fc Pugetl 11984 iDesert et all 11990( 1 . and loosely 
trace star formatio n rates on the scale of a galaxy 
([Roussel et al.|[200ll ). Spectral diagnostics have been de- 
veloped using the strengths of these featu res to distin- 
guish star formation f rom n uclear activity ([Genzel et al.l 
119981: lLaurent et all 12000( 1. PAH carriers are ap- 
parently ubiq uitous in high met a llicity star-formin; 
galaxies (e.g. iBoselli et al.l 119981 : iGenzel fc Cesarsk" 



120001: iDraine et all 12007( 1 . alt hough in low metallicity 
galaxies they may be absent ([Engelbracht et all 120051 : 
IDraine et al.l 12007). Calibration of physical quantities, 
such as star formation rate, from PAH bands is compli- 
cated by the relative importance of circumstcllar dust, 
and diffuse dust heated by the general radiation field 
or collisional impacts. These components may be sig- 
nificant for galaxies with low to mo derate star forma- 
tion rates and dust contents (see e.g. | Boselli et aDl2004 
iPerez-Gonzalez et alj|2006t ICalzetti et alj|2007[ ). 

To keep things simple, one can construct Mid-InfraRed 
(MIR) colours as opposed to calibrating star formation 
rates. The wealth of information in this spectral range 
prompted the use of MIR colours to trace star formation, 
dust temperature and nuclear activi ty (e.g. iDale et al.l 
120001 IBoselli et all 120031 : lFisherll2006l ). Normalizing the 
flux from a dust-sensitive band ( A rcs t > 6/im) by a 
stellar mass-sensitive band (tracing old stellar popula- 
tions, e.g. at 3.6/xm) produces a MIR colour highly 
sensitive to the level of dust emission per unit stellar 
mass. This is especially well correlated with morphol- 
ogy, tracing the warm dust emi ssion from star-forming 
spiral arms ([Pahre et al.l [2004a) . Indeed, in the local 
Universe, the MIR colou r can effectively predict mor- 
phology. iLi et al.1 ([2007H find 88% early-types below 
vL v (8)/z/L„(3.6) = 0.303 and 83% late-types above this 
division 1 . The combined sample exhibits a bi modal dis- 
tribution in MIR colour. iJohnson et al.l (|2007f ) also see 
this bimodality in MIR colours for galaxies in Hickson 
Compact Groups, and show that the fraction of MIR- 
passive galaxies is much higher in the more evolved (HI- 
poor and X-ray bright) groups. These results demon- 
strate the power of IRAC colours to separate truly pas- 
sive galaxies from star-forming or AGN-host galaxies 
without the dust-driven degeneracies of optical and UV 
studies. 

Our goal is to isolate truly passive galaxies and thus 
trace the build-up of the increasingly important popula- 
tion of passively evolving galaxies in the Universe, as a 
function of environment. In particular, we wish to study 
massive (M» > 2 X 10 10 M Q ) galaxies at z ~ 0.4, where 
the contributions of both passive and dusty star-forming 
galaxies are important, but can be confused with each 
other without due attention. Separating these popula- 
tions using the MIR colour will alleviate this problem. 
Strong evolution in the dust and star formation prop- 
erties of normal galaxies over the range < z < 0.7 is 
expected: At z > 0.7 spirals domina te the Luminous 
InfraRed Galaxy (LIRG) p o pulation dBell et ail 120051: 
iRowan-Robinson et all 120051 : Izheng et al J 12007( 1. whilst 
today galaxies of equivalent infrared luminosities are 
much rarer, and are typically interacting systems and 
mergers. Our ga laxy samp l e, sele cted from the CNOC2 
redshift survey (|Yee et al.1 12000( 1 . probes intermediate 
redshifts (z ~ 0.4). This redshift corresponds to a look- 
back time of ~ 3 — 4 Gyr, when the volume-averaged star 
formation rate in the Universe was significantly higher 
than it is tod ay ([Lilly et al.1 119961 : iMadau et ail 119981 : 
lHot)kinsl[200l . 

In Section[2]our sample is introduced, including archive 
IRAC data in our fields. From this stellar masses and 

1 Includes assumed stellar contribution to i/L„(8) of 
0.232^(3.6). 
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MIR colours are derived, a mass-selected sample is con- 
structed. Selection effects will be discussed in detail. 
Section [3] presents the MIR colours of z ~ 0.4 galaxies as 
a function of their stellar masses and morphologies. A 
tight Infrared Passive Sequence is identified in the MIR 
colour-mass plane as a tool to select passive galaxies. 
Section Q] examines the correlation of MIR colour with 
the EW[OII]A3727 emission line diagnostic, testing their 
relative sensitivity. In Section [5] we compute the fraction 
of massive galaxies with InfraRed Excess (IRE galaxies) 
as a function of stellar mass. The role of the group en- 
vironment is discussed in Section [6j wh ere we examine 
galax ies in our optical group catalogue (jCarlberg et al.1 
2001). Finally, our results are discussed in the context 
of our current understanding of galaxy evolution and the 
group environment in Section [7] In particular, by ap- 
plying a simple model to a mock galaxy catalogue, we 
focus on the possibility that downsizing is intimately re- 
lated to, if not entirely caused by, the growth of structure 
with time. Our conclusions are presented in Section [8l 
Throughout this paper a ACDM cosmology with Cl m = 
0.3, n A = 0.7 and H = lOO/i km s^ 1 Mpc^ 1 with h=0.75 
is assumed. 

2. SAMPLE 

The objective of this paper is to study MIR activity in 
galaxies as a function of their stellar mass and environ- 
ment. This is facilitated by the construction of a mass- 
selected sample in the redshift range 0.3 < z < 0.48, 
described in this section. 

2.1. CNOC2 Survey and Groups 

Our sample is based upon the CNOC2 redshift survey 
which is magnitude limited down to Rc ~ 23.2 in pho- 
tometry and totals ~1.5D° within four separate patches 
on the sky. Spectroscopic redshifts exis t for a large an d 
unbiased sample of Rc < 21.5 galaxies (|Yee et al.l l2000). 
Within the fields of 20 kinematica l ly sel ected groups 
from the catalogue of ICarlberg et al.l (|2001l ), the original 
UBVRcIc photometry and spectroscopy have been sup- 
plemented with targetted spec troscopy to Rc = 22 using 
the Magellan 6.5m telescope (|Wilman et al.ll2~005f) . and 
F775W-band HST ACS imaging. Galaxies with spectro- 
scopic redshifts z > 0.3 and ACS coverage have been 
visually morphologically classified by one of the authors 
(Augustus Oemler Jr, hereafter AO). 

The probability that a galaxy has a redshift is a 
function of its Rc-band magnitude, and is unbiased to 
Rc = 21.5. To study the statistical properties of galax- 
ies as a function of their stellar mass (and not Rc-band 
magnitude), we must statistically correct for this selec- 
tion function and magnitude limit. This is achieved by 
weighting galaxies by the inverse probability that they 
appear in the Rc < 21.5 spectroscopic sample. There 
are two components to this weight: 

• A weight Wr c to account for the fraction of galax- 
ies which are targetted spectroscopically (with 
measured redshifts) as a function of Rc-band mag- 
nitude. Since below Rc = 21.5 the success of 
redshift determination may be biased in favour of 
emission line galaxies, Wr c is set to for Rc > 
21.5. 



• A weight W z to account for the fraction of sample 
volume within which a galaxy would fall out of our 
Rc < 21.5 sample (effectively a V max correction). 
For a 0.3 < z < 0.48 sample, this means that for 
each Rc < 21.5 galaxy in the sample, W z = V(z = 
0.3 to z = 0.48)/V(z = 0.3 to z = z Um ) where V is 
the cosmological volume in our chosen cosmology, 
and zii m is the redshift at which the galaxy would 
have a magnitude Rc = 21.5. 

Then the final weight for each galaxy is simply: W = 
W Rc x W„. 

2.2. IRAC Photometry 

The Spitzer Space Telescope GTO program 64 (PI 
Fazio) included an IRAC (Infrared Array Camera, 
iFazio et alj 12004) scan-map survey of ~20'x30' regions 
in the centre of three of the four CNOC2 patches. The 
IRAC focal plane geometry translates each mapping into 
one mapping of bands 1 (3.6/im) and 3 (5.8/mi) and an- 
other mapping of bands 2 (4.5/im) and 4 (8.0/im), offset 
by 6.7'. 

The astrometrically and photometrically calibrated 
post-BCD (Post - Basic Calibrated Data) products were 
retrieved from the Spitzer archive. The pipeline ver- 
sion used was Sll.4.0. The processed images are clean, 
with no image artifacts, and provide stable source flux 
measurements and accurate sub-arcsecond astrometry. 
Apart from at the map edges, all sources receive ^500 s 
imaging time in all four bands (100 s frame-time at five 
dithered positions). The complete four band IRAC cov- 
erage has significant (> 50%) overlap with six ACS fields, 
where galaxies have also been morphologically classified. 
Shallow MIPS 24/im imaging is also available, but with 
lower sensitivity and poorer diffraction- limited resolution 
than the IRAC bands. 

For all sources detected in the 3.6/im frame, 
colours are measured using the SExtractor package 
(|Bertin fc Arnoutsl 119961 ) in double-image mode. They 
are derived from fluxes measured inside 3" radius aper- 
tures centred at the position of all 3.6/im detected ob- 
jects, in all IRAC bands. Detection requires at least 
four adjacent pixels which are at least 1.5-cr above the 
background level at 3.6/xm. Aperture corrections were 
computed to correct for light lost from the aperture due 
to the spatial extent of the PRF (Point Response Func- 
tion) 2 . Corrected measurements are used to compute 
flux ratios (colours) within the aperture. The SExtrac- 
tor FLUX_AUTO Kron-like measurement provides more 
stable total flux measurements, and thus is used to com- 
pute absolute galaxy luminosity at 3.6/im. We follow the 
practice of SWIRE 3 and compute the flux error from the 
measurement image itself, weighted by inverse variance 
using the coverage map. 

To match the 3.6/im selected catalogue to the CNOC2 
redshift catalogue, we apply the tmatch tool in iraf . 

2 Measured by combining larger aperture measurements for iso- 
lated stars with IRAC Data Handbook values. Corrections are 
1.184, 1.180, 1.249 and 1.404 for IRAC bands 1-4 respectively. 

3 swire.ipac.caltech.edu/swire/astronomers/publications 

4 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 
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None of the galaxies in our spectroscopic sample are lost 
from the IRAC sample due to lack of sensitivity except 
at the edge of the IRAC field. Visual inspection shows 
that CNOC2 galaxies with no match either lie outside 
the IRAC zone of full coverage, can be identified with 
sources improperly deblended by SExtractor, or lie in 
regions of particularly poor CNOC2 astrometry. The 
sample is then limited to those objects with unambiguous 
matches, with no apparent bias relative to the ~ 15% of 
CNOC2 galaxies with IRAC coverage which are lost at 
this stage. However, to maintain statistical accuracy the 
weight Wr c is modified to account for galaxies without 
spectroscopy or 3.6/im matches, as a function of Rq- 
band magnitude. At 3.6/im, the 2-cr depth is ~l/iJy (for 
a passive galaxy at z=0.4, this corresponds to a stellar 
mass of ~ 9 x 1O 8 M ). 

The goal is to construct a 0.3 < z < 0.48 mass-selected 
sample for which representative spectroscopy is available, 
and reaching required depth at 8/im to distinguish pas- 
sive galaxies from those with excess emission at 8/im (2-cr 
depth of ~10//Jy). Effectively, these criteria set the mass 
limit to 2 x 10 10 M Q , well above the 3.6/mi limit (§[13 
for more details). 

2.3. Deriving Stellar Masses and MIR Colours 
2.3.1. Stellar Mass 
Near InfraRed (NIR) light is a useful tracer of stel- 



mation history ( 


Aaronson et al.ll 19791 iRix & Riekdll993t 


iBell & de Jonsl 


2001). 3.6/mi flux has been converted 



into luminos i ty an d thence to stellar mass following 
iBalogh eT al. (2007). To give a brief overview of this pro- 
cedure: Observed-frame AB luminosities are computed 
for the assumed cosmology, and an empirical k-correction 
is applied to correct to a consistent rest-frame of a z = 0.4 
galaxy. In this paper, we will call this [3.6]. Secondary 
corrections are then applied to move into the K-band 
rest-frame. These corrections are calibrated for objects 
with both [3.6] and K-band photometry (not used in this 
paper). The mass to light ratio of each galaxy in K-band 
(M »/Lr) is estimated using a combination of two sim- 
ple | Bruzual &; Charlotl (|2003l ) models with a IChabrieil 
(|2003f) IMF. For blue galaxies (B - V < 0.4), a constant 
star formation model with dust extinction t v ~1 mag 
gives M*/L K = 0.2. For red galaxies (B - V > 1.0), a 
11.7-Gyr old, dust free single stellar population model 
gives M*/Lk = 0.7. For intermediate colours a linear 
interpolation betwee n these values i s used . For more de- 
tails, please refer to IBalogh et al.1 {2007). This means 
that the full range of galaxy types only exhibit a range 
in M*/Lk (and thus also M*/Lo.4[ 3 . 6 i) of 3.5. Therefore 
our results arc relatively insensitive to whether we se- 
lect in mass or 4 [3.6] luminosity. We prefer to select in 
mass as it is the more physical quantity and, thus, allows 
a comparison with the literature. 

2.3.2. MIR colour 

The two template spectra in the top panel of Figure [1] 
demonstrate the diversity and wealth of spectral informa- 
tion at MIR wavelengths. A passive, dust-free galaxy is 
likely to exhibit a spectrum similar to the SSP template 
(dashed line). With no dust emission, the infrared emis- 
sion originates solely from photospheres of the cold stel- 



lar population, primarily giant M-stars. Conversely, the 
spectrum of a dusty star forming galaxy should look more 
like that of M82 (solid line). Beyond A rcs t ~ 3 — 4/im, a 
typical dust emission spectrum with characteristic PAH 
features comes to dominate the M82 template. 

The lower panel of Figure [1] shows the transmission 
function of the IRAC bands at 3.6/im and 8/im, and 
the MIPS (Multi- Band Imaging Photometer for Spitzer, 
iRieke et al.ll2004T ) 24/im band. Each transmission func- 
tion is transformed to indicate the rest-frame wavelength 
domain probed for a galaxy located at z = 0.3,0.4 or 
0.5. It is clear that excess (enhanced beyond pure stel- 
lar photosphcric) emission can be revealed at both 8/im 
and 24/im for galaxies within the redshift range of in- 
terest. For 0.05 < z < 0.48 galaxies, the 8/im band 
captures the strong 6.2/im PAH feature, tracing star for- 
matio n (e.g. iPeeters et al.l 120041 iForster Schreiber et al.1 
2004). Crucially, stellar atmospheric emission drops off 
significantly from 8/im to 24/im: Thus detection of pas- 
sive galaxies, and thence accurate separation of passive 
from star-forming or AGN-host galaxies, requires much 
shallower exposures at 8/im than 24/im. A colour con- 
structed from 8/im and 3.6/im has the leverage on the 
galaxy spectrum to make this distinction. In this paper 
we focus on IRAC data for these reasons. 

For the sake of simplicity, we compute k-corrected MIR 
colours of our galaxies using a linear combination of the 
tw o simple templates i n FigurefT] (similar to the approach 
of iHuang et al.ll2007h . We define our colour to be the 
flux-ratio [f(8)/f(3.6)] 5 . For each galaxy, the quantity 
77 is measured: This is the required fractional contribu- 
tion from the M82 template to match the observed-frame 
z [f(8)/ f(3.6)] 6 . Our k-corrections are based on the as- 
sumption that r] (approximately dust to stellar flux ra- 
tio), would be the same regardless of the redshift of the 
galaxy (and thus the part of the spectrum sampled), 77 
is computed as: 



n 



z [f(8)/f(3.6)] - z [f(8)/f(3.6)](SSP) 
? [f(8)/f(3.6)](M82) - z [f(8)/f(3.6)](SSP) 



(1) 



To minimize the k-corrections z [f(8)/ f(3.6)] is trans- 
formed to the equivalent colour for a z=0.4 galaxy, and 
limit ourselves to galaxies in the narrow redshift range 
0.3 < z < 0.48, an epoch > 3Gyr ago. The k-corrected 
[f(8)/f(3.6)] colour is computed as: 

a4 [f(8)/f(3.6)] = 77X a4 [f(8)/f(3.6)](M82)+(l-r ? )x °' 4 [f(8)/f(3.6)](e 

(2) 

A pure SSP k-correction (< 0.07 z [f(8)/ f(3.6)]) in- 
creases by < 10% going from ages of 8 Gyr to 1 Gyr, 
and 10 — 20% going from solar to 3x solar metallicity. 
For 7/ = 1 (M82-like), 4 [f(8)/ f(3.6)]=2.17 (~ 7 times 
the stellar photospheric continuum emission within the 
IRAC 8/im band). Whilst k-corrections of high 77 galax- 
ies are more uncertain due to the assumed M82 template, 
k-corrections remain relatively small within our redshift 
range (< 50% for 77 = 1). In this paper we will only sepa- 

5 for a more traditional colour, simply compute [3.6]-[8] = 
+2.51ogio([f(8)/f(3.6)]) 

6 The prefix supersci pt z is used to indicate the band-pass 
(Blanton & Rowcis 2007); i.e. the redshift of a galaxy for which 
the k-correction is defined to be zero. 
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Fig. 1. — Our two MIR spectral templates are illustrated in the top panel. A 4Gyr old Simple Stellar Population traces pure stellar 
atmospheric emiss ion from old sta rs (SSP, dashed line). This is generated using PEGASE2 (Fioc & Rocca-Volmcrange 1997), with solar 
metallicity, and a Salpeter (1955J) IM F with upper mass limit of 120 Mq. The other template is the ISO-SWS spectrum of the dusty 
star forming galaxy M82 (solid line, Forstcr Schrcibcr ct al. 2001). The data is limited to A > 2.4/im, and so we extrapolate down to 
1.9/jm by combining it with the spectrum obtained from a composite stellar population model (constant star formation rate, evolving 
metallicity, 1 Gyr since start of star formation, Salpeter IMF with upper mass limit of 120 Mq). The two templates are normalized to 
match at A = 2.2/im. Beyond ~ 3fim, the spectrum of M82 is dominated by dust emission in the form of PAH features (< 12/^m) and a 
warm dust emission continuum (> 12/xm). The bottom panel shows the effective transmission functions of Spitzer IRAC bands at 3.6/^m 
and 8fim, and the MIPS 24/xm band in the rest-frames of galaxies at z = (solid), and z = 0.3 (long-dashed), 0.4 (short-dashed) and 0.5 
(dotted), corresponding to the redshift range of our galaxies. Whilst 3.6/im traces stellar emission, 8(im and 24/^m trace PAH and warm 
dust emission respectively at z ~ 0.4. 



rate dusty galaxies emitting excess light at infrared wave- 
lengths from passive galaxies, and will not draw strong 
conclusions from precise MIR colours of dusty galaxies 
which are sensitive to the applied k-correction. To this 
extent, we believe that we are justified in applying these 
rough k-corrections. 

2.4. A Mass-Selected Sample 

In this paper, our analysis is limited to the subsample 
of galaxies with known redshifts 0.3 < z < 0.48, and with 
full IRAC coverage. Of these, the 333 galaxies with stel- 
lar masses M* > 2 x 10 10 M Q are selected, corresponding 
to a luminosity of a4 [3.6]=-20.89 for passively evolving 
galaxies, or a flux of f(3.6^m) > 20/iJy, well above the 
detection limit. 

This stellar mass limit is actually defined by the depth 
of the 8/xm data, since we require this colour informa- 
tion for our analysis. In our mass-limited sample, 313 
of the 333 galaxies are detected (2a) at 8/im, and the 
remainder have > 1.8er upper limits that are sufficient 
to identify them as passively evolving galaxies, charac- 



terised by colours 0A [f(8)/ f(3.6)] < 0.5 (see § [SI for a 
justification of this definition). 

The requirement that all galaxies have a spectroscopic 
redshift introduces sensitivity to the optical (Rc) selec- 
tion function. To account for selection bias, the spectro- 
scopic selection weight W is applied ( §2.1|) . This weight 
includes the factor W z which applies a volume correction 
to correct for galaxies which fall out of the Rc < 21.5 
sample at z Vlm < 0.48. Within aM,>2x 1O 1O M sam- 
ple, the lowest mass and highest redshift galaxies with 
high mass to light ratios will not make this magnitude 
cut. Nonetheless, the weighted sample is still representa- 
tive: A passively evolving M* — 2 x 10 10 M Q galaxy can 
be sampled with Rc < 21.5 up to a redshift of zi; m ~ 0.38 
(~ 36% of the volume sampled) and receives a weight 
W z = 1./0.36 = 2.78. There do exist four R c < 21.5 
galaxies with unusually high mass to optical light ra- 
tios. These can more easily be lost from the sample 
due to their extremely red optical- NIR colours. Such red 
colours can exist where there is extreme dust extinction, 
or a population of AGB stars boosting NIR luminosity. 



() 



However, they are too few to significantly influence our 
statistics. Our sample of 333 galaxies includes eleven 
M» >2x 1O 1O M galaxies with R c > 21.5. These are 
assigned zero weight, but it is worth noting that the IRE 
fraction of these optically faint galaxies is comparable to 
that of other galaxies. In practice, our results are not 
strongly sensitive to the statistical weights. 

3. BIMODAL MIR COLOURS AND A PASSIVE SEQUENCE 

In this section the distribution of the MIR colour 
4 [f(8) / f(3.6)], and its dependence on galaxy stellar mass 
M* is investigated. Figure [2] shows how the 333 galaxies 
in our M« > 2x 10 10 M Q , 0.3 < z < 0.48 sample populate 
the °- 4 [f(8)/f(3.6)]-M* "colour-mass" plane, along with 
the 25 th , 50 th and 75 th percentiles of the distribution. 

The distribution shifts to lower 04 [f(8)/ f(3.6)] with 
higher mass. Dust flux per unit stellar mass can be ap- 
proximated by 4 [f(8)/ f(3.6)] 50 % - 0.308 where the stel- 
lar atmospheric contribution to 4 [f(8)/ f(3.6)] is esti- 
mated by the SSP template to be 4 [f(8)/ f(3.6)](SSP) = 
0.308 (dotted line). The median relation between 
4 [f(8)/ f(3.6)] - 0.308 and stellar mass can be well fit by 
a linear relation in log space: log\o{ °' 4 [f(8)/ f(3.6)] 50 % — 
0.308) = 8.8 ± 0.2 - (0.85 ± 0.02)logi ( M*/ M Q ). This 
relation is significantly shallower than would be ex- 
pected for no relation between dust flux and stellar mass 
(°- 4 [f(8)/f(3.6)]-°- 4 [f(8)/f(3.6)](SSP) ~ ^). However, 

this is sensitive to the exact value of 4 [f (8)/ f (3.6)] (SSP) 
chosen. 

We also examine whether the strong correlation be- 
tween MIR a ctivity and morpho l ogy observed i n the lo- 
cal Universe (jPahre et al.ll2004al : ILi et al.|[2007l) extends 
to our intermediate redshift sample. 31 of the galaxies 
in our sample are within the ACS fields and have been 
morphologically classified by AO. These are indicated in 
Figure [2j 

Morphologically classified early-type (Elliptical and 
SO) and late-type (mainly spiral) galaxies exhibit very 
different MIR colour distributions, as seen locally by 
ILi et alj (|2007f ). In fact, the overall distribution in 
4 [f(8)/ f(3.6)] is bimodal in nature, with the two peaks 
in MIR colour almost equivalent to these two morpho- 
logical classes. This bimodality has bee n obseryed in 
other (local) optical- NIR selected samples (|Li et al.ll2007l : 
iJohnson et alJ l2007h whilst it is less obvious in shal- 
low 8/im-selected s amples where early -types are more 
difficult to detect (jHuang et all I2007D . The bimodal- 
ity of CNOC2 galaxies is demonstrated in Figure 02 
split into two bins of stellar mass. Within each mass 
range there exists a significant population of galaxies 
with 04 [f(8)/f(3.6)] < 0.5 and then a gap before the 
onset of a second population of galaxies emitting more 
strongly at 8/zm. Whilst the full distribution seems 
to shift to higher 04 [f(8)/ f(3.6)] at lower mass, the 
value 04 [f(8)/ f(3.6)] = 0.5 is within the gap across 
our mass range. Thus we adopt this value to divide 
the sample: Everything with 4 [f(8)/ f(3.6)] > 0.5 we 
call an InfraRed-Excess (IRE) galaxy. This division 
at 0.5 excludes galaxies from the passive population if 
they have 44% more 8/Ltm flux than the early-type lo- 
cus. Contribution from circumstellar dust (jPiovan et alJ 
2003 1 iTemi et al.ll2008ft or silicate emission from mass- 
loss (|Bressan et al.ll2007j) are unlikely to produce such a 
strong boost to the 8/im flux. 



The dash-filled histogram in Figure [3] indicates the 
04 [f(8)/ f(3.6)] distribution for morphologically early- 
type galaxies. The tight distribution of high mass 
early-types is mirrored in the narrow locus of almost 
constant colour in the colour-mass plane (Figure [2]). 
The 11 early-types (ignoring the single upper limit) 
with M* > 6 x 1O 1O M have mean colour and scat- 
ter 4 [f(8)/f(3.6)] ~ 0.347 ±0.03 (horizontal dashed line 
in Figure E|) ■ The scatter is consistent with the photo- 
metric errors. We call this the Infrared Passive Se- 
quence (IPS), which is not far offset from the tem- 
plate SSP colour ( 04 [f(8)/f(3.6)] = 0.308, horizontal 
dotted line). The peak value and scatter of the IPS 
for the full population of M* > 6 x 10 10 Mq galaxies 
is 4 [f(8)/f(3.6)] - 0.34 ±0.05. The tight scatter of 
the IPS is reminis cent of the optical red sequence (e.g. 
iBower et al.l[l9 92h Indeed the tightness of the IPS is 
due to the emission from stellar atmospheres of cold 
stars in passive galaxies, which is relatively little sen- 
sitivity to age or metallicity. Mid Infrared wavelengths 
are highly sensitive to dust emission; thus the Infrared 
Passive Sequence does not suffer from the degeneracy be- 
tween dusty red galaxies and passive red galaxies as by 
a red sequence defined at optical wavelengths. 

The second, broader peak at 4 [f(8)/ f(3.6)] > 0.5 con- 
tains IRE galaxies with excess (non-stellar) emission at 
8/im. This excess emission traces some kind of activ- 
ity, i.e. star formation or AGN. Figures and show 
that all known late-type galaxies (dot-filled histogram in 
Figure [3|) in the sample inhabit this peak. Two morpho- 
logically classified early-type galaxies also exhibit excess 
emission at 8/im. Excess MIR emission in some early- 
type galaxies is known to exist, and correla tes with low 
levels of star formation or A GN activity fe.g. lPahre et al.l 
l2004lAlBres~san et al.ll200l . 

4. THE 0,4 [f(8) / f(3.6)]-EW[OII] PLANE 

It is useful to compare the MIR colour 4 [f (8)/ f (3.6)] 
with an independent, well-studied indicator of activ- 
ity (primarily star formation), the equivalent width of 
the [OII]A3727A emission line. [Oil] is most commonly 
used at z>0.3 where the HaA6563A emission line is red- 
shifted out of the optical window. By normalizing the 
emission flux by the continuum level, the sensitivity to 
flux calibration and dust absorption is reduced, produc- 
ing a normalized quantity more easily comparable to 
°- 4 [f(8)/f(3.6)]. 

Figure S] shows the location of our M, >2x 10 10 M Q 
and 0.3 < z < 0.48 galaxies (excluding 14 galaxies where 
EW[OII] could not be measured) in the 4 [f (8)/ f (3.6)]- 
EW[OII] plane, including 25 th , 50 th and 75 th percentiles 
and morphologically classified galaxies overplotted with 
their errors. 

There is a strong correlation between EW[OII] and 
4 [f(8)/f(3.6)]. Both diagnostics trace galaxy activity, 
and primarily star formation, and so a strong correla- 
tion is expected. A value of EW[OII]= 2.1A (dotted 
vertical line) corresponds to the median galaxy with 
04 [f(8)/f(3.6)] = 0.5 (our empirical division in MIR 
colour, solid horizontal line), which is below the [Oil] 
detection limit of most spectroscopic surveys. Indeed, 
we are only able to identify this value for a large statisti- 
cal sample: The median error on EW[OII] for individual 
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Fig. 2.— The °' 4 [f(8)/ f(3.6)]- M* "colour-mass" plane for our sample of 333 M, > 2 X 1O 1O M , 0.3 < z < 0.48 galaxies. For the few 
galaxies without 2<r 8/^m detections, 2a upper limits on a4 [f(8)/ f(3.6)j are indicated. 25 th , 50 th and 75 th percentiles of the 4 [f (8)/ f (3.6)] 
distribution are overplotted (red, solid and dashed lines). These are unweighted, and computed for running bins of 50 galaxies in M*, and 
undetected galaxies are considered to be located at their measured values as opposed to upper limit values. Larger symbols are used for 
galaxies in the joint IRAC-ACS sample for which galaxies are morphologically classified (see key). For morphologically classified galaxies, 
measurement errors on 0,4 [f(8)/ f(3.6)] are shown, which are typical of the overall population. Early-type galaxies cluster along a tight 
sequence at 4 [f(8)/f(3.6)] ~ 0.35, consistent with pure stellar photospheric emission (horizontal dashed line). We call this the Infrared 
Passive Sequence (IPS), and define InfraRed Excess (IRE) galaxies as those with 0,4 [f(8)/ f(3.6)] > 0.5 (horizontal solid line). The SSP 
template is indicated by the horizontal dotted line ( 0,4 [f(8)/ f(3.6)] = 0.308). 



gala xies in our sample is 2 . 4 A . 

In IBalogh et all (|2007f ). we divided our sample into 
passive and [OII]-strong galaxies at EW[OII] = 10A. This 
is the equivalent of a division at [f(8)/f(3.6)] = 
1.31 for the median galaxy, selecting only the tail of 
galaxies with the most ongo ing activity. A div ision 
at EW[OII]=5A (as used in iWilman etafl l2005h cor- 
responds to 4 [f(8)/f(3.6)] = 0.92. Therefore with a 
reasonable photometric depth at A rcst ^6/^m, it is pos- 
sible to pick up levels of activity not possible without 
extremely deep optical spectroscopy. This inevitably 
leads to a higher fraction of galaxies defined as non- 
passive. Many spiral galaxies (mainly of type Sa-Sbc) 
with 4 [f(8)/f(3.6)] > 0.5 have no detectable [Oil] emis- 
sion. Based upon a simple cut in EW[OIIl, these would 
be classified as pas sive spirals (e.g. lPoggianti et aT1ll999l : 
IBalogh et al.1 20021. It is now clear that most such galax- 
ies are not truly passive (in the sense of the IPS galaxies): 
instead there is evidence at 8/im of some ongoing activ- 
ity. This indicates that the morphological transforma- 
tion is more intimately linked to (and concurrent with) 
the complete cessation of star formation than previously 
believed. 

The distribution of galaxies in the 04 [f(8)/ f(3.6)]- 



EW[OII] plane exhibits a lot of scatter around the me- 
dian relation. Both EW[OII] and a4 [f(8)/ f(3.6)] are 
highly sensitive to the presence of dust (through ex- 
tinction of the [Oil] emission by dust, and abundance 
of PAH carriers), which is itself de pendent on the 
galax y mass and star formation rate dGiovane li et al.l 
[19951: iWang fc Heckmanl fl996t iMasters et all 120031 : 
iBrinchmann et aLll2004h . Thus whilst MIR diagnostics 
are better at tracing activity in high mass galaxies, ac- 
tivity in lower mass galaxies may be better traced with 
optical diagnostics (especially at significantly sub-solar 
metallicity and thus very low mass where PAH feature s 
disappear, lEn gclbrach t et aI1l2005l : [Draine et al.l l2007). 
However, within the mass range traced here, the median 
relation for EW[OII] goes as ~ M~ - 8 , similar to the de- 
pendence of 4 [f(8)/ f(3.6)] (once corrected for the stel- 
lar atmospheric contribution). To compute EW[OII], the 
[Oil] line flux is normalized by the continuum at A3727A, 
which will lead to a shallower slope than would a normal- 
izion by mass. It is likely that this effect counters any 
steeper relation expected due to different dust dependen- 
cies. Therefore, other physical parameters are necessary 
to drive the scatter in the 4 [f(8)/ f(3.6)]-EW[OII] plane, 
together with large measurement errors. 
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Fig. 3.— The weighted distribution of the MIR colour °- 4 [f(8)/f(3.6)] for 0.3 < z < 0.48 galaxies with M, >6x 10 10 M (top) and 
M» = 2 - 6 X 10 10 M Q (bottom), binned in lo g-space and with poissonian error-bars. The few galaxies without 2<r detections at 8fim are 
excluded. Since these are all constrained at the > 1.8cr level to lie at °- 4 [f(8)/ f(3.6)] < 0.5, their inclusion would not destroy the bimodality. 
Overplotted are the distributions for morphologically classified early type (dash-filled) and late type (dot-filled) galaxies, renormalized to 
give the same maximum value. Our criteria for an InfraRed Excess galaxy ( 4 [f(8)/ f(3.6)] = 0.5) is indicated with the vertical dashed 
line, and effectively separates the peak of the Infrared Passive Sequence (IPS) ( 0,4 [f(8)/ f(3.6)] < 0.5) from the InfraRed Excess (IRE) 
galaxy population. 

a4 [f(8)/f(3.6)] down to very 
star formation and nuclear 



There exists a strange population of galaxies extend- 
ing to high EW[OII] at low 4 [f(8)/f(3.6)]. These galax- 
ies diverge from the usual situation in massive galaxies: 
Namely, that infrared data traces activity to much lower 
levels than emission lines. For two of the four galax- 
ies with 4 [f (8)/ f (3.6)] < 0.5 and EW[OII]> 20A, the 
rcdshift depends upon a single emission line and might 
therefore be incorrect. Another possibility is that the 
[Oil] emission originates in a low mass/metallicity galaxy 
which is falsely matched to a bright IRAC source. Visual 
checks suggest this is unlikely, but cannot be certain. Al- 
ternatively, these galaxies might be unusually bright at 
3.6/im for their mass or metallicity. This may occur in 
metal-poor high mass galaxies (weird) or galaxies with a 
strong and very hot AGN component contributing most 
significantly to the 3.6Mm-ban d (and AGN activ ity can 
also boost the [Oil] emission, lYan et al.l ()2006l )). The 
only morphologically classified galaxy at high EW[OII] 
and low 0,4 [f (8)/ f (3. 6)] is a merger galaxy. This could be 
explained by either of the above scenarios, or a third sce- 
nario in which a merger-driven starburst produces such 
a hard radiation field that it is capable of destroying the 
PAH carriers. Whatever their origin, these galaxies form 
a puzzling, possibly interesting, but small minority. 

5. FREQUENCY OF INFRARED-EXCESS MASSIVE 
GALAXIES 

Within our representative sample of 0.3 < z < 0.48 
and M, > 2x 10 10 M Q galaxies, we now measure the 
fraction of InfraRed Excess galaxies, f(IRE). 



Activity is traced with 
low levels at which both 
activity play a role. We will not distinguish between 
these two forms of activity in this paper, instead merely 
computing the total fraction of IRE galaxies. How- 
ever, any galaxy which has gas accreting onto a super- 
massive bl ack hole is likely to be simultaneously form- 
ing stars (iBrinchmann et all l2004t iBerta et all 120071 : 
iShi et alJ 120071 ). Hence such galaxies are correctly dis- 
tinguished from the passive population. 

To measure f(IRE), the fraction of galaxies with 
4 [f (8)/ f(3.6)] > 0.5 is computed, with galaxies weighted 
to account for the magnitude-dependent selection func- 
tion (Section I2.4|) . The application of this weight does 
not significantly change the results. To place errors on 
f(IRE), we perform bootstrap resampling of the sample. 
This is done by creating 1000 independent samples of the 
same size as the true sample: Within each sample, values 
of °- 4 [f(8)/ f(3.6)] are taken randomly and often repeat- 
edly from the true sample, and scattered using a gaussian 
kernel of width equal to the error on 4 [f (8)/ f(3.6)]. The 
standard deviation of f(IRE) from the bootstrap samples 
is taken for the error on f(IRE). 

Figure E](a) (solid points and line) shows f(IRE) 
binned by stellar mass for the global population. Bins 
are not equally populated. Instead they are defined such 
that there is an approximately equal number of group 
galaxies per bin (§E|. 4 [f(8)/ f(3.6)] declines with M*, 
an effect seen previously in the median and percentiles. 
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Fig. 4.— The °- 4 [f(8)/f(3.6)] - EW[OII] plane for 319 M* > 2 X 10 10 M Q and 0.3 < z < 0.48 galaxies. 25"\ 50" 1 and 75 th percentiles 
of the °- 4 [f(8)/ f(3.6)] distribution are overplotted, computed over running bins of 50 galaxies in EW[OII] (red, solid and dashed lines). 
Undetected galaxies (at the la level in [Oil] or the 2a level at 8^tm) are plotted at their la EW[OII] and 2a °' 4 [f(8)/ f(3.6)] upper 
limit levels (arrows), whilst the percentiles are computed based upon measured [Oil] and 8/^m flux levels: on a statistical basis the 50*^ 
percentile should therefore be representative of the true median relation (both measured [Oil] and 8/im flux can be negative for non- 
detections). Our empirical division (° [f(8)/f(3.6)] = 0.5) is overplotted (solid horizontal line). There is a strong correlation between 
EW[OII] and °- 4 [f(8)/ f(3.6)], inferring that for the average galaxy both diagnostics trace the same phenomenon (primarily star formation). 
The median galaxy at our empirical division 4 [f(8)/ f(3.6)] = 0.5 has EW[OII]= 2.1 A (dotted vertical line). Larger symbols are used for 
morphologically classified galaxies as in Figure [2] with measurement errors typical of the overall population. 



However this is only evident in the highest mass bin. 

To test for effects of cosmic variance, this exercise is 
repeated independently for each of the three CNOC2 
patches with IRAC coverage 7 . The patch to patch off- 
sets in f(IRE) are < 3%, compatible with the bootstrap 
errors, and indicating that cosmic variance is not a dom- 
inant source of error. Therefore this result is robust. 

The crossover mass (M cr , where f(IRE) = 0.5: dotted 
horizontal line), is M cr ~ 1.3 x 10 M Q (dashed ver- 
tical line). iHopkins et alJ ((2007) compiled estimates of 
M cr (their "M tr ") at different redshift. These are scaled 
to a Chabrier IMF and our cosmological parameters for 
comparison with our measurement. Defining M cr as the 
stellar mass at which 50% of galaxies are optically blue or 
have a high specific star formation rate (calibrated using 
emission lines), they measure M cr ~ 1.6 — 3.8 x 10 10 M Q 
at z ~ 0.4, almost an order of magnitude lower than our 
estimate. This probably reflects our sensitivity to low 
levels of star formation in high mass galaxies. However, 
a definition of 50% morphologically late-type galaxies re- 
sults in an estimate of M cr - 3 - 8 x 10 10 M at z ~ 0.4, 

7 CNOC2 consists of four patches spaced at intervals of ~ 90° 
in RA, of which three arc partially mapped by IRAC. 



much more compatible with our value within errors. This 
suggests that morphological selection is identifying more 
or less the same population as infrared selection. This is 
interesting, as it tells us that morphological transforma- 
tion usually goes hand in hand with the final truncation 
of star formation in massive galaxies. 

6. THE GROUP ENVIRONMENT 

We will now examine how the fraction of InfraRed Ex- 
cess galaxies depends upon the group environment by 
computing f(IRE) for group galaxies as a function of stel- 
lar mass. Groups were detected as overdensities in red- 
shift space by applying a fr iends-of-friends algor ithm to 
the CNOC2 redshift survey ( |Carlberg et al.ll2001h . Mem- 
bership has been reassigned after inclusion of additional 
spectroscopy an d application of the membership algo- 
rithm defined in lWilman et alJ (|2005h 8 . To ensure most 
group galaxies are not interlopers (i.e. nearby galax- 
ies not bound to the group at the time of observation) 
the sample is limited to those within 500 kpc of the 
luminosity-weighted group centre. 500 kpc corresponds 
to the virial radius for a ~ 360 km s _1 group at z=QA. 

8 The group aspect ratio, b, is set to 5. 
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Fig. 5. — (a) The dependence of f(IRE) (the fraction of InfraRed Excess galaxies, defined by 4 [f(8)/ f(3.6)] > 0.5) for galaxies at 
0.3 < z < 0.48, divided into four bins of stellar mass (M,). This is computed for the full (global) sample (solid line, solid points) 
and the group galaxy sample (dashed line, open points). Points are plotted at the median M, of galaxies in each bin, and errors on 
f(IRE) are computed using a bootstrapping method (see text). There is a significant deficit of InfraRed Excess galaxies in groups, which 
nonetheless mirror the stellar mass dependence of the global population (f(IRE) decreases at high mass). A "crossover mass" (M cr ) at 
which f(IRE) = 0.5 (dotted horizontal line) is defined: For the global sample M C r ~ 1.3 X 10 11 Mq (dashed vertical line), (b) As (a), 
except f(IRE) is plotted against 0,4 [3.6] luminosity for galaxies in a luminosity-limited sample. 



Our decision to select group galaxies within a fixed ra- 
dial aperture instead of a scaled aperture (e.g. the virial 
radius) is made because measurements of velocity disper- 
sions are uncertain, and known t o be heavily biased for 
syste ms with < 30 members (e.g. IZabludoff k, Mulchaevi 
1998). 78/333 galaxies from the mass-selected sample 
match these group criteria. Of these 59% are in groups 
with measured velocity dispersions, Oi n t r < 360 km s _1 ; 
80% are in groups with <7i n t r < 500 km s _1 ; and 6% are 
in the cluster g226 (ointr ~ 850 km s _1 ). 

The dashed line and open symbols in Figure 0(a) show 
how f(IRE) depends upon stellar mass for group galax- 
ies only. Each bin contains 19 or 20 group galaxies by 
design, and appears to show a deficit of IRE galaxies 
with respect to the global population. This is not highly 



significant in individual mass bins, but is consistent at 
all masses. Thus we estimate the combined significance 
for the whole sample. In groups, the combined value of 
f(IRE) is 0.447. As the mass distribution of group galax- 
ies is biased to high mass relative to the global sample, 
the global sample is resampled to match the group sam- 
ple on mass. For each group galaxy, a random galaxy 
is selected from the nearest 10 galaxies in stellar mass 
(from the global sample, other than itself, and with no re- 
peats). This resampling process is repeated 10000 times, 
each time computing f(IRE) for the matched sample. On 
only 56 out of 10000 occasions is f(IRE) < 0.447. This 
corresponds to a 0.56% likelihood (~ 2.5cr from a one- 
tailed gaussian) that the low value of f(IRE) in groups 
arises by chance. 
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The lower value of f(IRE) in groups independent of the 
stellar mass indicates that activity in M* > 2 x 10 10 M 
galaxies is often suppressed by the group environment. 
As these groups are predominantly low velocity disper- 
sion, often unvirialised systems, this means that massive 
galaxies are already affected by the low mass halo en- 
vironment. The mechanism of suppression seems to be 
particularly active in the highest stellar mass bin ( M* > 
1.27 x 10 11 Mq) where in groups f(IRE) = 0.12 ± 0.07. 
In this mass bin the group galaxies also dominate the 
passive global population: Out of a total of 43 galaxies 
(25 passive), 20 are in our group sample (17 passive). 
Hence the most massive galaxies seem to become passive 
in the group environment. This result is even more re- 
markable when it is considered that the group sample is 
incomplete due to incompleteness effects in the CNOC2 
redshift survey, and so the global sample should include 
more group galaxies than have been identified. 

It is prudent to repeat this exercise as a function of 
4 [3.6] luminosity to ensure this result does not some- 
how depend upon our mass calibration. Figure 0(b) 
shows how f(IRE) depends upon 4 [3.6] luminosity for a 
luminosity-selected sample (°' 4 [3.6]< —20.89). Luminos- 
ity selection allows an additional 180 M* < 2 x 10 10 M 
galaxies to make it into the sample with lower than max- 
imum mass to light ratios. Effectively this means that 
the overall value of f(IRE) is higher for the luminosity- 
selected sample (0.735±0.018) than for the mass-selected 
sample (0.621 ± 0.030). The trends mirror those seen as 
a function of mass with the exception that in the bin of 
highest luminosity the global f(IRE) is still quite high 
(0.645±0.047). This is because IRE galaxies are more 
likely to be found in the high luminosity bin than in the 
high mass bin, due to their lower mass to light ratios. In 
groups the dominance of passive galaxies is so strong as 
to mitigate this effect. 

7. DISCUSSION 

MIR colour is a highly sensitive tracer of activity in 
galaxies. This has allowed us to separate truly pas- 
sive galaxies (typically morphologically early-type) from 
galaxies with InfraRed Excess at 8/im (IRE galaxies). In 
this way we avoid confusion between galaxies with low 
but still ongoing star formation and those which have 
shutdown their star formation altogether (probably co- 
incident with the morphological transformation). Other 
studies examine the evolution in fractio ns of highly star 
forming galaxies (LIRGs, ULIRGs etc: iBell et all 120051: 
Le Floc'h et al.l 120051: [Hammer et al.l 120051: ICaputi et al.1 
2006tlDaddi et al.ll2007f) . These fractions decline steeply 
with cosmic time, with contributions from generally de- 
clining star formation rates in most galaxies, and com- 
plete suppression in some. 

Our main result is that a strong and highly signifi- 
cant deficit in the fraction of InfraRed Excess galaxies, 
f(IRE), exists in optically selected, mostly low velocity 
dispersion groups at z~0.4, and is independent of stellar 
mass. The suppression of infrared activity in lo w redshift 
Hickson Compact groups (|Johnson et al.l[2007h therefore 
extends to the "loose group" regime at z~0.4. 

Our group data indicate a dramatic drop in f(IRE) 
at stellar masses M* > 1O 11 M0. This coincides with 
the stellar mass at which a galaxy's own halo usually ex- 
ceeds Mhaio ~ 5 x 10 12 Mq, which is also the mass of 



the Local Group (|Li fc White! l2007f) and typical of the 
lower mass groups in our sample (McGee et al., submit- 
ted). Therefore galaxies with M* > 10 11 M are typi- 
cally "central" galaxies in groups, which is why many of 
our most massive galaxies are members of known groups. 
The low mass group environment suppresses star forma- 
tion in both central massive galaxies and their satellites. 
The strongest effect is in the central galaxies. 

In this paper, we do not fully constrain the mecha- 
nism^) responsible for suppressing activity in groups. 
Fo r a discussion of the menagerie of possibilities we refer 
to iMoran et al.l (|2007l . and references therein). The abil- 
ity to detect low levels of activity has important implica- 
tions for the interpretation of pass ive spirals as identified 
purely from their [ OH] emission dPoggianti et al.l [1999; 
iBalogh et ail 120021 : IMoran et all \2QM). Our strongest 
direct constraint on the mechanism(s) at work is that it 
(they) must be active in low mass group haloes. This 
probably excludes processes such as ram-pressure strip- 
ping which prefer much higher dens ity environments 
ijGunn fe G"ot^ll972tlQuilis et al.ll2000D . 

Relating the suppression of star formation to a halo 
mass threshold is theoretically attractive, as the physics 
of gas accretion and feedback are more likely to re- 
late to the gravitationally bound material within the 
halo than to the stellar mass of the galaxy itself (e.g. 
iDekel fe Birnboimll2006T ). For example, AGN feedback 
processes might be especially efficient in galaxies located 
inside group-sized haloes. Brightest Group Galaxies 
(BGGs) more frequently contain radio-loud AGN than 
similarly massive non-BGGs. This is likely to be related 
to the cooling of the hot Intra-Group Medi um (IGM) 
onto the BGG, fuelling a radio-loud AGN (jBest et al.l 
[2001 . 

Our data alone is not sufficient to examine how the 
suppression of star formation within group-mass haloes 
contributes to the global population in a statistical man- 
ner. Incompleteness of the CNOC2 redshift survey 
and uncertainties in the construction of the group cat- 
alogue and galaxy membership inevitably leads to in- 
completeness of our group galaxy sample in a way which 
may depend upon galaxy stellar mass. To better esti- 
mate membership on a statistical basis, we use a sim- 
ulated galaxy catalogue in which stellar mass and em- 
bedding halo mass can be directly evaluated for each 
galaxy at z=0A. We choos e the "mock" galaxy cat- 
alogue of Bo wer et al.l (|2006l ) which is based upon the 
Millenium Simulatio n, a simulated ACDM Universe with 
500/ifoo Mpc sides (jSpringel et all l2005fh These dark 
matter haloes are populated with mock galaxies using 
a semi-analytic technique. Only the positions (environ- 
ment) and stellar masses of these galaxies are used. The 
mock galaxies are assumed to be in groups where they are 
located within haloes Mhaio > 5 x 1O 12 M . Mock group 
galaxies are then assigned a probability that they are 
IRE galaxies consistent with the data: i.e. f(IRE) = 0.5 
for the M* < 10 11 M galaxies and f(IRE) = 0.1 for 
the M* > 1O 11 M galaxies (as seen in Figure 03(a)). 
Combined with the rest of the mock galaxy population 
("field" galaxies), the dependence of f(IRE) upon stellar 
mass is compared to the CNOC2 global trend. By set- 
ting the fraction of IRE field galaxies to 0.9, regardless 
of their stellar mass, the global trend is well reproduced. 
This is demonstrated in Figure [6] (a) in which the ob- 
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served global trend for CNOC2 galaxies (solid points) is 
extremely well matched by the model (solid line). That 
the fraction of field galaxies which need to be passive is 
low (0.1) and independent of stellar mass lends support 
to our hypothesis that the halo mass is more important 
than the stellar mass for the suppression of star forma- 
tion in galaxies. 

Figure[6](a) also shows how f(IRE) depends upon stel- 
lar mass for the mock galaxy population evaluated at 
z = 0, 0.7, and 1. Remarkably, the model exhibits a 
downsizing behaviour, such that the crossover mass M cr 
decreases at lower redshifts. This is shown in Figure [5] 
(b) with the CNOC2 datapoint at z=0A. The only rea- 
son for this evolution is that the fraction of galaxies in 
groups increases, as a function of stellar mass, from z = 1 
to z = 0. Suppression of galaxy activity in groups, plus 
the gravitational growth of structure in the Universe, can 
explain downsizing trends in M cr . 

Also in Figure [5] (b) the crossover ma s ses co mputed in 
different redshift bins by iBundv et all {2006) are over- 
plotted (reproduced from their Figure 9: M cr is called 
the ^transition mass" in that paper). These measure- 
ments are based on the DEE P2 redshift survey, an d con- 
tribute to the compilation of Hopkins ct al. (2007). M cr 
is computed for 50% late- type morphologies, 50% galax- 
ies blue in (U-B) colour, and 50% galaxies with [OII]- 
inferred SFR above O.2M0yr _1 . That more morpholog- 
ically classified late-type galaxies are detected, resulting 
in a higher value of M cr , is fully in line with our esti- 
mate at z=0.4 using a sensitive MIR colour, well cor- 
related with morphology. Whilst our simple model can 
qualitati vely reproduc e a do wnsizing scenario, compari- 
son with Bun dv et al.l (|2006h suggests that the redshift 
dependence is probably too shallow. With better con- 
straints at different redshift, a model in which star forma- 
tion is not immediately suppressed upon incorporation of 
a galaxy into a massive halo might serve to steepen this 
evolution. 

8. CONCLUSIONS 

From the CNOC2 redshift survey, we have constructed 
a mass-selected sample of galaxies with known spec- 
troscopic redshifts 0.3 < z < 0.48, and stellar masses 
M, >2x 1O 1O M . Whilst the sample is not fully com- 
plete, the selection is well understood, and we correct for 
incompleteness using a simple weighting scheme. Within 
this sample, the fraction of galaxies with excess emis- 
sion at 8/im is computed, which traces star formation or 
nuclear activity. For galaxies in this sample, the Mid In- 
fraRed (MIR) colour 4 [f(8)/ f(3.6)] has been measured. 
This is the flux ratio of the 8/im to 3.6/im IRAC bands, 
k-corrected to a consistent rest-frame of a z=0.4 galaxy. 
Based upon this sample we have shown: 

• Infrared Passive Sequence: Old stellar popula- 
tions (such as found in passively evolving galaxies 
in which star formation has been suppressed) ex- 
hibit very predictable colours in the ~ 2 — 7/xm 
spectral region (rest-frame), tracing the stellar at- 
mospheres of giant M-stars. This produces a 
tight Infrared Passive Sequence (IPS) in the 
04 [f(8)/f(3.6)] colour. The IPS does not suffer 
from the optical red sequence degeneracy between 
dusty red galaxies and passive red galaxies. The 



mean IPS colour ( 4 [f(8)/ f(3.6)] = 0.34 at M» > 
6 x 1O 1O M0) is consistent with an old stellar pop- 
ulation, and the tight scatter (0.03 — 0.05 at high 
mass) is consistent with photometric errors. This 
leaves little room for large variations in metallicity, 
or additional contributions at ~ 6/im from circum- 
stellar or diffuse dust. 

• MIR Colour Bimodality: The 4 [f(8)/ f(3.6)] 
distribution is bimodal, consisting of the IPS and a 
population of galaxies with significant excess emis- 
sion at 8/im (InfraRed-Excess Galaxies). These two 
populations can be divided at 04 [f(8)/ f(3.6)] = 
0.5. 

• Morphologies and infrared emission: 15 out 

of 17 morphologically classified early- type (ellip- 
tical or SO) galaxies have °' 4 [f(8)/ f(3.6)] colours 
consistent with the IPS, whilst no classified late- 
type galaxy has a4 [f(8)/ f(3.6)] < 0.5. This is 
consistent with known correlations at low red- 
shift between morphology and M IR colours (e.g. 
iPahre et aT1l2004aHLi et al.ll2007[ K 

• Tracing low level activity in the MIR: There 
is a strong correlation between 04 [f(8)/ f(3.6)] 
and EW[OII]. However, our division at 
°' 4 [f(8)/ f(3.6)] = 0.5 corresponds to a division at 
EW[OII]= 2.LA for the median galaxy. This is 
below the [Oil] detection limit of most intermedi- 
ate to high redshift surveys. This shows that MIR 
diagnostics trace low levels of activity at z ~ 0.4. 
Consequently, dividing at 04 [f(8)/ f(3.6)] = 0.5 
selects a high fraction of InfraRed-Excess (IRE) 
galaxies, f(IRE), and the crossover mass at which 
f(IRE) = 0.5, M cr ~ 1.3 x 10 n M o , is also high 
with respect to optical studies. 

• Consequences for selection of passive spirals: 

There exist spiral galaxies in our sample with un- 
detected [Oil] emission but significant MIR excess 
emission (°- 4 [f(8)/f(3.6)] > 0.5). Thus, selection 
of passive spi rals based upon a simple EW[OII] 
divis i on (e.g. Poggian ti et alJ Il999t iBalogh et al.l 
l2002t iMoran et al.ll2006l ) will suffer contamination 
by dusty, low level star formation, mostly early- 
type spirals. 

• Strong suppression of activity for massive 
galaxies in low mass groups: In the z ~ 0.4 

optically-selected group environment there is a 
strong and robust (2.5u) deficit in f(IRE) com- 
pared with global values, independent of stellar 
mass. This is true across our range in stellar mass 
M* > 2 x 10 10 M Q , but the suppression appears to 
be particularly strong at high mass M* > 10 11 M Q , 
where the fraction of active group galaxies falls to 
f(IRE) = 0.12 ±0.07. 

• Suppression in groups can drive the stel- 
lar mass dependence, and structure growth 
can drive downsizing: The global dependence of 
f(IRE) on stellar mass at z=0.4 can be accurately 
reproduced by a simple model in which suppression 
mainly occurs in groups, and is particularly strong 
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Fig. 6. — (a) The dependence of f(IRE) on stellar mass M, for the global population of 0.3 < z < 0.48 CNOC2 galaxies compared 
with our simple model predictions at Z. = 0.4. The model predictions are also computed at z = 0, 0.7, and 1. Galaxies in haloes with 
Mhalo — ^ x 10 12 M© are considered to be in groups. The stellar mass dependence can be reproduced by a model in which suppression 
of galaxy activity is driven primarily by the group environment, (b) Evolution of the crossover mass M cr (where 50% of all galaxies are 
active, computed in various ways). Our data point ( M cr at z=0A computed using f(IRE) = 0.5 in CNOC2 ) is compared with the model 
predictions at z = 0,0.4,0.7, and 1, and DEEP2 crossover masses computed in different redshift bins by|Bundy ct al. (2006, B06) for 
50% late-type morphologies, 50% galaxies blue in (U-B) colour, and 50% [Oil] -inferred SFR above 0.2Mq yr" 1 . Our simple model can 
qualitatively reproduce a downsizing scenario, although the redshift dependence is shallower than observed by B06. 



for M* > 10 11 M group galaxies. Furthermore, 
structure growth in the Universe drives more galax- 
ies into group-sized haloes with cosmic time. This 
automatically leads to a downsizing phenomenon in 
the global population, in which the crossover mass 
M cr evolves to lower masses at lower redshift. 
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